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Abstract 
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We make use of the superoperator technique. 
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I. INTRODUCTION 



Master equations, or evolution equations for reduced density operators are of fundamental 
importance in the treatment of open systems in quantum theory. They are particularly 
relevant in the field of quantum optics [I|,0] , where generally one has interest in studying the 
evolution of confined subsystems (a single mode of the electromagnetic field, for instance) 
"coupled" to external reservoirs. This may model a number of (phase-sensitive or not) 
processes, such as the decay of a single mode of the field confined in a lossy cavity JE], as 
well as amplification processes HH. Interesting applications of phase- sensitive amplification, 
for example, are the reduction of noise in lasers through the injection of squeezed vacuum 
Q, as well as noise- free amplification via the two-photon correlated-emission-laser ||. 

Nevertheless, complete operator solutions of master equations are not usually presented. 
Only simple cases such as the decay of a single mode of the field into a vacuum state 
reservoir (T = 0) or stationary regime cases are the ones normally treated H0. Instead, 
master equations are normally transformed into c-number equations either in the number 
or coherent state representations, for instance. It would be therefore interesting to obtain a 
direct solution of the master equation while in its operator form. 

II. THE METHOD 

In this paper we show how to solve exactly master equations of the form 




where 

Lx(a,a))p = 2apa) — a)ap — pa) a (2) 
L 2 (a, a))p = 2a) pa — aa) p — pad) (3) 
L 3 (a,a^)p = 2apa — a 2 p — pa 2 , (4) 

and 
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1/4 (a, a))p = 2a) pa) — a) 2 p — pa) 2 , (5) 

being a and a) the bosonic annihilation and creation operators. The 7s are in general 

complex parameters which may represent gain and decay. For the density operator to be 
Hermitean, it is necessary that 

7s = |7s|e> = 74- ( 6 ) 
We may now apply the unitary transformation 

p = S(t)p'&(t) (7) 
where is the squeeze operator 

m ^ (£* - <*!) , (8) 

with £ = rexp(iip). If we now substitute Eq.(|7|) into Eq. ([]]), we obtain an equation for the 
transformed density operator p', or 

I^XjTM^' (9) 

where 

6 = S0aS(£) = pa + va), (10) 

with 

p = cosh(r), v = — sinh(r) exp(iip). (11) 

We proceed in rewriting Eq.(|S|) in terms of the original creation and annihilation opera- 
tors as 

do' 

= Li (a, a))p' (71/i 2 + 72M 2 + HTsle^ + pu*\j 3 \e-^) + 
L 2 (a, a))p' ( 7l |i/| 2 + l2 p 2 + HTsle^ + A*^* I Ta I e - * ) + 
L 3 (a, ( 7lyU z/* + j 2 pu* + p 2 \-f 3 \e^ + {y*f\ lz \e^) + 
L 4 (a, a))pt Upv + j 2 pu + ^73^ + p 2 ^^) . (12) 



We note that one can choose <p an d r in such a way that the last two terms in Eq.(|12| 
become zero. The appropriate choices for these parameters are 



<P = -<, 



(13) 



and 



tanh(2r) 



2|7a| 



7i + 72 

With the above choice of parameters, then we can write Eq. (|i~2"D as 



dl_ 
dt 



7iLi(a, a f ) +7 2 L 2 (a, a 1 ) 



(14) 



(15) 



where 



7i = 7i/" 2 + 72^| 2 + /H73|e J<?i + //i/* 1 73 1 e 



(16) 



and 



72 = 7iM 2 + 72/i 2 + ^l73|e l<?i + //i/*|7 3 |e" 



(17) 



It follows from the relations above that 



7i - 72 = 7i - 72, 



(18) 



and 



7i + 72 = (7i + 72) sech(2r). 



(19) 



We now write the new master equation fll5|) in a more convenient form [[7] 

?t = + j 2 + J 3 - 27 2 ) p', 
where we have defined the following super-operators 



(20) 



(21) 



J 2 f> = 2^ 2 a?p'a 



(22) 



and 



JzP = —(71 + 7 2 )(a t ap' + p'a'a) 



(23) 



From the equations above, we can write the formal solution of Eq.([20|) as 



P 



'(*) = exp(-27 2 t) exp + J 2 + J 3 ) t] p'(0). 



(24) 



It is not difficult to show that the superoperators Ji, J 2 and J3 obey the following 
commutation relations IHI 



J2, J\ 



P 



47i72 



7i + 72 



— J 3 - 47i7 2 



P- 



(25) 



4-y p' = -2 (71 + 72) Ap', 



(26) 



and 



J21 J3 



p' = 2 (7! + 72 ) J 2 p'. 



(27) 



In order to disentangle the exponential in Eq . (|24|) , we propose the ansatz 



p'(t) = exp(-2 72 t) exp[/ 3 (t)] exp[f 2 (t)J 2 ] exp[f Q (t)J 3 ] exp\fi(t)Ji(t)]pf(p) 



(28) 



By inserting p'(t) in Eq.(p8|) above into equation in Eq. (p0[) , we obtain the following 
system of differential equations for the functions ff. 



dfo 47172 

-77- + ~ — —J2 - 

at 71 + 7 2 



(29) 



dt 



iexp[2(7i + 72)] = 1, 



df2 n df j, _ v . „ 2 1 

Tjt" 7 ^ 47172 ^ 2 = lj 



(30) 
(31) 



and 



^ - 47172/2 = 0. (32) 

Although the system is a non-linear one, its solution is rather straightforward. In order 
to have the condition p'(t = 0) = p'(0) satisfied, we should have /i(0) = (i — 0, 1, 2, 3) as 
initial conditions for the set of Eqs.(p9|)-(|3~2|). The result is 



k 1 {K + T(t)\ , s 

fo = — —t + — — In K -L), 33 

7i + 72 7i+72 V K J 



* = * = iA < 34 » 



and 



h = 272t - In (^^) , (35) 

where T(t) = (1 - exp[-2«t])7 2 . 

By finally defining the dimensionless superoperators 

L-p = hpb , L + p = b pb, , L 3 p = Pbp + pPb + p, (36) 

where the superoperators L_, L + and L 3 obey the commutation relations [L_,L + ]p = L 3 p 
and [L3, L±]p = ±2L±p, we can write the solution of ([!]) in the form, 

La 



p{t) = e Kt e«+ r w 



Kt pK +r(t)- L + 



« + r(t) 



5l r(i) L 



e^^w^-p(O), (37) 



where in (36) we have defined 

b = S(i)aS\i) = pa-va\ (38) 

which should be compared with Eq. (|T0|). Note that if in Eq. ( pTf ) we set the parameters 
£ = and 72 = we recover the usual solution for a dissipative cavity at zero temperature 
as well as for the phase insensitive case. 

We have therefore obtained the full solution of the master equation for p' [see Eq. (|20|)1, 
given by (|37|). We remind that we still have to apply the transformation in Eq.(|7]) in order 
to recover the original density operator p. 
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It would be appropriate now to show which interesting cases could be easily treated 
by employing our solution. For instance, if we make 71 = 7 (n + l)/2, 72 = 7I/2, 73 = 
-7M72 and 74 = — 7M/2, we will obtain the standard equation describing the decay of 
a bosonic field into a phase-sensitive reservoir 0. If on the other hand 71 = An/2 and 
72 = A (n+ l)/2, but keeping 73 = —AM*/2 and 74 = -iM/2, we will have the standard 
equation describing phase- sensitive amplification The constants 7 and A represent decay 
and gain, respectively, and M in) are connected with phase-sensitive (or not) reservoir 
fluctuations. Each of this two cases are relevant in the problems of field decay and in the 
reduction of noise in lasers [|J. We remark that our operator solution turns unnecessary the 
convertion of the master equation into sometimes cumbersome c-number equations. 

Let us finally consider as an example an initial density matrix of the form 

p(0) = 5(O|0)(0|£t(O, (39) 

or a squeezed vacuum state. After inserting it into fl3"?D it is easy to show that the density 
operator of the field at a time t may be written as 

r(t) 



\m,0(m,£\, (40) 



where = S(^)\m) are the squeezed number states P,P^|. We may now to follow the 



state's evolution in phase space, using for instance the Q-function [0, defined as 



Q(t) = -((3\p(t)\(3), (41) 

7T 



where \/3) is a coherent state. In our case the Q-function is given by 

r(t) 



K + T(t) 



m,0\ 2 , (42) 



with 



KO = expH/3172) ]T G nm (-/ry\ (43) 

n=0 



and where 



G, 



(sSb) 1/2 (!^)" +m/2 El .„ ^.^l,,), if «, m even 

(-l)" + "' /2 - 3/2 (=W) 1/2 (^) , " + ™ ,/2 - 1 E« p , + ,„ ( ^,{^ ) , if »,m odd 
otherwise 

(44) 

In figure 1 it is illustrated the Q-function in ([Sf) for different times. At t = we have 
the Q function of the initial squeezed vacuum state. As time goes on, thermal fluctuations 
of the reservoir cause a "spread" of the Q function, associated to the increase in the quadra- 
ture noise. Nevertheless, due to the phase-sensitive properties of the reservoir, the noise 
assymmetry characteristic of squeezed states is somehow preserved, as it may be seen in 
figure 1. 

III. CONCLUSIONS 

To summarize, we have presented an alternative way of treating problems involving 
certain types of master equations. Our operator solution may be useful to retrieve any sort 
of information related to the evolution of sub-systems having specific initial conditions. This 
contrasts with the usual approaches, where only partial information, such as mean values of 
amplitudes and/or diffusion coefficients is normally obtained. 

Note added in proof: While preparing the answer to the referee, we became aware of 
the paper by Dung and Knoll [11] where they also solve Eq. (jl]) but using Fokker-Planck 



equations instead. We would like to stress that our method of solving the master equation 
(d) is considerably more straightforward and simple. 
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FIGURES 






FIG. 1. Q function of the cavity field initially in a squeezed vacuum state, a) at nt = 0, b) at 
Kt = 0.5, and c) at nt = 1.0. It has been taken £ = r = 0.7 real, and 72 = 1. 
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